Phase separation driven by a fluctuating two-dimensional self-affine potential field 
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We study phase separation in a system of hard-core particles driven by a fluctuating two- 
dimensional self-affine potential landscape which evolves through Kardar-Parisi- Zhang (KPZ) dy- 
namics. We find that particles tend to cluster together on a length scale which grows in time. The 
final phase-separated steady state is characterized by an unusual cusp singularity in the scaled cor- 
relation function and a broad distribution for the order parameter. Unlike the one-dimensional case 
studied earlier, the cluster-size distribution is asymmetric between particles and holes, reflecting the 
broken reflection symmetry of the KPZ dynamics, and has a contribution from an infinite cluster 
in addition to a power law part. A study of the surface in terms of coarse-grained depth variables 
helps understand many of these features. 



I. INTRODUCTION 

The behaviour of a scalar field driven by a fluctuat- 
ing force field depends strongly on the correlations of the 
driving field in space and time. A well known example 
is the passive scalar problem in fiuid mechanics, where 
one asks for the behavior of a passive field as it is ad- 
vected by a turbulent fluid flow [Q. Even in situations 
where the driving force field has simpler correlations, the 
passive scalar field can show interesting, and sometimes 
unexpected, behaviour. In particular, while in the fiuid 
context an initial local concentration of passive particles 
typically spreads out in space, in other types of situations 
an initially randomly distributed set of particles may be 
driven into a state with large-scale clustering. In this 
paper, we study one such example. 

We consider a force field which is derived from a fluc- 
tuating potential, and ask for its effect on a system of 
particles which do not interact with each other except 
through hard-core exclusion. The problem is then tan- 
tamount to the dynamics of hard-core particles which 
reside on a fluctuating surface and are driven downwards 
by gravity along local slopes. An especially interesting 
case arises when the surface is self-affine with a power-law 
divergence in the height correlation function as a func- 
tion of the separation. In such cases, surface roughening 
strongly affects the clustering of particles and can lead 
to new types of states [||,^. We study a stochastically 
evolving two-dimensional surface governed by Kardar- 
Parisi-Zhang (KPZ) dynamics [Q and show that surface 
fluctuations bring about large-scale clustering of parti- 
cles, akin to phase separation. The particles are taken to 



be random walkers with an excluded volume constraint, 
diffusing in the dynamic potential landscape defined by 
the local height h{r,t) at a base position r = {rx,ry) at 
time t. In addition, the particles are assumed to be suf- 
ficiently massive that the dynamics operates effectively 
at zero temperature; this means that particles only move 
locally downwards at a fixed rate, subject to the con- 
ditions that the local slope is favourable and the target 
site is unoccupied. This model is a generalization to two 
dimensions of the 1-d model studied in ^ where it was 
found that the particles reach a phase separated state 
with unusual properties arising from strong surface fiuc- 
tuations. Since the nature and strength of fiuctuations 
depends strongly on the dimensionality, it is important 
to see to what extent these features survive in higher di- 
mensions. This is one of the principal aims of this paper. 

Let us summarize the main results. Our numerical 
simulations support the argument that the particle den- 
sity exhibits phase ordering over a characteristic time- 
dependent coarsening length scale C{t) ^ t^^^ set by sur- 
face fluctuations Here z is the dynamical exponent 
for surface fluctuations; for the KPZ surface with d ^ 2, 
it is known that z ~ 1.6 The pair-correlation func- 
tion is found to depend on the separation scaled by C(t), 
as is characteristic of a system undergoing phase order- 
ing, but the scaling function has a cusp near the origin, 
unlike usual phase separating systems j^]. The cluster 
size distribution in the steady state has a power law de- 
cay N(s) ~ s^'^ for s <^ LF' , where the exponent r is 
different for particles and holes. In addition, unlike the 
1-d case, there is a distinct contribution to N{s) from 
an 'infinite' cluster which contains a finite fraction of the 
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total number of particles. An understanding of these re- 
sults can be gained from a study of the surface itself. 
To this end, it proves useful to distinguish between re- 
gions where the height is less than or more than a fixed 
reference level. This is incorporated in a Coarse-grained 
Depth (CD) model of the surface HQ, which itself ex- 
hibits phase separation. Many of its properties reflect 
the underlying asymmetry of the KPZ surface. Finally, 
in both the shding particle and depth models, the order 
parameter has a distribution that remains broad in the 
thermodynamic limit. 

We introduce the model in Section II. In Section III, 
we present results for two-point correlation functions, 
cluster-size distributions and the order-parameter distri- 
bution, while Section IV is the conclusion. 

II. MODEL 

We study an autonomously evolving two-dimensional 
surface, whose behaviour over large length and time 
scales is described by the Kardar-Parisi-Zhang equation 

i 

— ^iyy^h + XiVhf+r,{v,t) (1) 

with short range correlations for the noise, 
{r]{r,t)r]{r',t')) cx S'^{r ~ r')S{t - t'). The KPZ sur- 
face \\\ d — 2 is known to be self-afRne in steady state: 
((/i(0) — h{r))'^) ^ r^^ where x is the roughness expo- 
nent. For the KPZ problem with d = 2, x - 0.4 [|,|). 
The nonlinear term in Eq.(l) breaks h — > ~h symmetry, 
and this has important consequences for the models we 
study. 

We simulate the surface through a discrete solid-on- 
solid (SOS) algorithm [H, where the height difference be- 
tween nearest neighbour (NN) points on a square lattice 
is maintained at ±1. A point is selected at random and 
its height is increased by 2 units with probability p+ if all 
four of its NN points are at greater height, and decreased 
with probability p_ if all four are at a lower height. Oth- 
erwise, the site is not updated. It is believed that the 
asymptotic properties of this single step model are the 
same as those of the (2-|-l)-dimensional KPZ equation, 
though this has not been proved. In our simulations we 
chose p+ = and p_ = 1, so that we have an 'evaporat- 
ing', rather than a growing surface whose average height 
decreases in time. 

The Sliding Particle (SP) model is defined as follows. 
Particles are initially distributed at random on surface 
sites with no more than one particle per site, the overall 
particle density being p. The external force field which 
drives the particles acts downwards, so that particles 
tend to move in the same direction as the average surface 
height, while holes tend to move upwards, opposite to the 
direction of surface motion. In a microstep, a randomly 



chosen particle attempts a move to a randomly chosen 
neighbouring site. If the local slope is favourable and the 
target site is unoccupied, the move is made, otherwise the 
particle stays at the same site. One MC step is counted 
when, on average, all surface sites and all particles have 
been updated in a random sequence. We verified that a 
faster rate of updating for particles with respect to sur- 
face updates does not change the qualitative results. We 
did our simulations on square lattices with p = \ and 
sizes ranging from L = 32 to L = 256. We used peri- 
odic boundary conditions for both the surface and the 
particles which reside on it. 

It is helpful to introduce two sets of discrete Ising- 
like spin variables to characterise the particle and surface 
configurations. For the particle configuration, we define 
cr(r) — 2n(r) — 1 where n(r) is the local occupation in- 
dex, i.e., n(r) = 1 if there is a particle at r and zero oth- 
erwise. Further, to characterise the height fluctuations 
of the surface configuration, we define a Coarse-grained 
Depth model in which we categorize sites according to 
whether they are above or below a certain fixed height. 
To this end, we define s(r) = —sgn{h{r) — ho) where ho 
is a chosen reference level. We call the set {cr(r)} as SP 
spins and the set {s(r)} as CD spins. In order to have a 
rough correspondence with the half-filled case p = ^ oi 
the SP model, we take ho to be a spatial average over the 
configuration of surface heights at that instant. 

We expect a correlation between SP and CD spins as 
local slopes in the surface guide the particles towards the 
local minima, so that over sufficiently large time scales, 
particles are expected to preferentially cluster in low- 
depth regions with predominantly positive CD spins. If 
particles were to occupy the lowest available positions so 
as to minimise the total potential energy, we would have 
a close match of the CT(r) and s(r). However, this 'ground 
state' is never reached as the surface reconfigures itself 
before particle rearrangements can occur. 



III. RESULTS 

A. Correlation Functions 

The SP-CD correspondence gives an insight into the 
extent of clustering of the particles. For example, a hill 
(a region with negative s(r)) of linear extension ^ is ex- 
pected to overturn in typical time t ~ causing clus- 
tering of particles as they fall into the valleys that form. 
Since the largest hill that overturns within time t has 
typical size C{t) ~ t^^^, this defines the characteristic 
length scale for particle clustering This is verified 
in a numerical study of the pair spin correlation func- 
tion C{r,t) = {a{0,t)a{r,t)) in the SP model. As we 
see in Fig.l, this function can be fitted into a scaling 
form C(r,t) = f{r/C{t)) with C{t) ~ t^/^ and z = 1.6, 
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which is close to the value of the dynamical exponent for 
the (2+1) dimensional KPZ surface ||^. This dynamic 
scaling form is characteristic of a phase ordering system. 
However, the scaling function has a cusp singularity at 
small argument, unlike usual coarsening systems where 
the function is linear for small arguments pj. 
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FIG. 1. Equal-time two-point spin correlation for parti- 
cles sliding down a 2-d KPZ surface of linear size 256. There 
is a scaling collapse when plotted as a function of the scaled 
variable r/t'^ with (j> — 0.6, indicating coarsening. Each point 
represents an average over 100 different time histories. 

In steady state, the system size L replaces C{t) as the 
relevant length scale, and a similar scaling form is ex- 
pected for the correlation function. To quantify the de- 
viation from linearity at small arguments, we studied the 
structure factor in the steady state, 

S{\^,L) = jd'^rC{r,L)e-'^\ (2) 

As shown in Fig. 2, the direction-averaged struc- 
ture factor is well described by the form S'(fc, L) — 
I - co + L'^g{kL), where cq ~ 0.4 and giq) - g~(2+") 
at large q. The first term arises from a short-distance 
analytic contribution (1 — co)d{r) which adds on to the 
scaling part /(r/L) of C(r,L). We find a ~ 0.38(4), 
a pronounced difference from typical coarsening systems 
where a — 1 (Porod Law) ]To[ |. The non-Porod form im- 
plies a cusp in the real space scaling function at small 
argument, 

f{x)cico-cix°' + ... (3) 

with X — r / L. 

We have also studied the two-point correlation func- 
tion in the CD model, and find a similar behaviour, with 
a cusp exponent a ~ 0.43(3) In the next subsection we 
will show that within the independent interval approxi- 
mation, the cusp exponent a of the CD model is equal 
to the roughness exponent x- 
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FIG. 2. Steady state pair correlation in the SP model 
for various system sizes L. The finite size scaling form 
C{r,L) = f{r/L) illustrated by the data collapse is is char- 
acteristic of phase separation. Inset: The scaled structure 
factor g{q) plotted against the scaled wave vector q = kL 
with L = 128 shows that g{q) ~ with a ~ 0.38(4), 

implying a cusp singularity in the real space correlation func- 
tion. 

B. Cluster Size Distributions 

Let us define a cluster as a set of like SP spins (parti- 
cles or holes), each of which is a nearest neighbor of at 
least one other like spin in that cluster. A study of the 
size distribution of clusters shows (Fig. 3) that the num- 
ber of connected clusters N±{s) with s particles (holes) 
has a power-law decay for s L^: 

N±{s) L^s-^^ ; s<L2 

where r+ ~ 2.2 and r_ ~ 2.0. This power-law distri- 
bution of cluster sizes, reminiscent of critical systems, is 
another characteristic feature of the unusual phase sepa- 
rated state under study — it occurs in one dimension as 
well There are, however, two important differences 
from the 1-d case. First, there is a marked difference in 
the powers r+ and t_ for particle (ctj = -1-1) and hole 
(ctj — —1) clusters, which reflects the asymmetry of the 
Si = +1 and Sj = — 1 clusters in the CD model as dis- 
cussed below. Second, for the particle cluster distribution 
we find that there is an additional large s contribution 
of the form foois/L^ — yo) peaked at yo ~ 0.4. For finite 
system sizes, foo is somewhat broad, but with increasing 
L, the width of the peak narrows down. We find that 
the area under the peak is unity, implying that there is a 
single very large cluster in every configuration. Also, the 
fraction of particles contained in this cluster is 2yo ~ 0.8. 

The power-law distribution of cluster sizes in the SP 
model has its counterpart in the CD model too. It is 
known that when a rough surface is intersected by a hor- 
izontal plane, the distribution of areas enclosed by the 
closed contours of intersection has the power-law form 
7V(s) ~ s-^' for s < L^, where r* = 2 - | |l|^. A 
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typical contour encloses several other contours, and the 
overall structure is scale invariant [^. For the KPZ sur- 
face one should moreover distinguish between contours 
whose inner perimeter sites lie above or below the cut, 
as there is no a priori reason for both the distributions 
to be identical. In fact, numerical simulations show that 
cluster size distributions for CD spins (which are closely 
related to area distributions), indeed has an asymmetry 
between positive and negative spins (Fig. 4). Both the 
distributions follow power-law decays, but with different 
exponents. The exponent for negative CD spins is nu- 
merically close to r*, while that for positive spins is sig- 
nificantly larger. Morover, there is a contribution from 
an infinite cluster of positive spins, as for the particle 
distribution in the SP problem. 
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FIG. 3. Size distribution of cluster sizes for -I- spins (par- 
ticles) and — spins (holes) for the SP model. The data shows 
a power-law decay with exponent ~ 2.2 and a well-separated 
'infinite' cluster for particles, and a lower value of the expo- 
nent ~ 2.0 for holes. We used L — 128 and averaged over 100 
histories. 

Within the CD model, we argue that power laws in 
cluster distributions of a different sort are related to the 
occurrence of a cusp in the scaled correlation function. 
For a self-afline surface with roughness exponent it 
is known that the probability P{1) that the surface first 
returns to its starting height /i(x = 0) = ho after mov- 
ing a distance I along an arbitrary linear direction has 
a power-law decay at small I: P{1) ^ ;~(2-x) for / <C £ 
|7|, p]Jl4|] . Each such segment defines a linear cluster of 
CD spins Si of the same sign, along the linear cut. Now 
let us make the independent interval approximation (IIA) 
in which the lengths of such segments are taken to be in- 
dependent random variables. This enables the Laplace 
transforms of the cluster size distribution and pair cor- 
relation function for CD spins (defined relative to ho) to 
be related along any linear cut. The nonstandard feature 
is that the mean cluster size diverges as L ^ cxd, but the 
correlation function can be calculated as in 0. It has 



the scaling form 

C*(r,L)~l-a(y)x + .. ; r<i {IIA) (5) 

Comparing with Eq.||, we see that the IIA predicts that 
the cusp exponent a in the CD model is equal to the 
roughness exponent %. The KPZ value x is quite close to 
the measured value of a for the CD model, and also to 
that for the SP model. The behaviours of both models 
on large scales of distance and time appear to be similar 
in 2-d, even though the microscopic configurations of the 
two match only roughly. To get an idea of the latter, we 
monitored the overlap index O — ((T(r)s(r)), and found 
that O varies between ~ 0.38 when the frequency of up- 
dates of particles and surface is equal, and a saturation 
value ~ 0.6 as the ratio between the two is increased. 
However, there was no corresponding significant change 
in the numerical values of the exponents a and r. 
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FIG. 4. Domain size distribution for -I- (valley) and — (hill) 
spins for the CD model. The exponents for the power-law 
parts have values ~ 2.2 and 1.85 respectively. There is a 
distinct contribution to the distribution of -I- spins from an 
infinite cluster. We used L — 128 and averaged over 100 
histories. 



C. Order Parameter 

For a system with conserved magnetisation like the SP 
model, an appropriate quantity to characterise the or- 
dered state is the steady state average of the magnitude 
of the Fourier components of the density Q, defined as 
g(k) = (|L-2X;r"(r)e"*-''|) where k = ^{n^,ny) and 
k| < TT. Taking the magnitude guarantees that (3(k) re- 
ceives the same contribution from all configurations that 
can be obtained from each other by translational shifts 
jisf . In Fig. 5, we plot Q(k) along the (1,0) direction 
for four different lattice sizes. The sequence of curves 
suggests that for any fixed, finite k, Q{k = 27m/L,0) 
approaches zero as L increases. However, if n is held 
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fixed so that |k| approaches zero as L — > cxo, the corre- 
sponding Q approaches a finite Hmit. Of this set of Q's, 
the largest is Q* = Q{^,0), and provides the simplest 
characterization of the order. 

I 1 1 1 1 1 1 1 

L = 32 + 
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FIG. 6. Probability distribution of the order parameter 
Q* for three lattice sizes for the SP problem. 

Fluctuations in the SP model are unusually strong and 
are reflected in the probability distribution of Q* (Fig. 
6). The distribution seems to remain broad and approach 
a limit as L increases, with (Q*) ~ 0.07 and variance 
~ 0.02. Similar studies of the corresponding CD model 
shows that there is a broad distribution there as well, and 
we find (Q*) ~ 0.11. The fact that the RMS fluctuation 
of the order parameter does not vanish in the large-size 
limit is an unusual and characteristic feature of the or- 
dered state of this model. Note that fluctuations which 
drive Q* towards zero need not take the system into a 
disordered state. A study of the temporal fluctuations 
of Q* in one dimension showed that a low value of Q* 



occurs simultaneosly with an increase of the value of Q 
for another close-by n, so that the state retains macro- 
scopic order, but with a few more coexisting macroscopic 
domains . A more complete characterization of the or- 
der in two dimensions would thus involve finding the joint 
probability V[Q(^^, ^^)] as i ^ oo, but this has not 
been attempted here. 

It is interesting to compare the distribution P{Q*) for 
this model with the corresponding quantity for a more 
familiar system such as the 2-d Ising model evolving un- 
der conserved (Kawasaki) dynamics. For the ferromag- 
netic Ising model with equal numbers of up and down 
spins, the ordered state in a finite system consists of 
strips of width ^ \L parallel to either x oi y axes, with 
equal probability. Thus the distribution P(Q*) (say, with 
Q* = <3(t"' ^) ) ^^^^ would have two peaks, one 
at a non-zero value of Q* corresponding to strips forming 
along the x-direction, and another peak at Q* = which 
corresponds to strips forming in the y-direction. In our 
case, however, the peak at the origin is absent as the in- 
terface between phases is much more diffuse and there is 
no strip formation. In this sense, the ordering observed 
in this model is very different from that in traditional 
equilibrium models. 

IV. CONCLUSION 

To conclude, we have studied a model where a fluctu- 
ating self-affine surface drives a set of downward-drifting 
particles to a phase separated state. The phase separa- 
tion has several unusual characteristics arising from the 
presence of strong fluctuations which survive even in the 
thermodynamic limit. These include a power-law distri- 
bution of cluster sizes, a cusp in the scaled pair corre- 
lation function and a finite width of the order param- 
eter distribution. These features mirror the properties 
of a coarse-grained height model of surface fiuctuations, 
where they follow directly from the self-affine nature of 
surface fluctuations. In particular, the asymmetry in the 
cluster size distribution for particles and holes is a con- 
sequence of the fact that KPZ growth breaks up-down 
symmetry. It would be interesting to see if a similar 
asymmetry exists in other related quantities such as the 
distribution of lengths of contours and the areas enclosed 
by them. 

It would also be interesting to study the sliding parti- 
cle problem on a fluctuating Edwards- Wilkinson surface, 
where particle-hole symmetry should be restored in the 
half-filled case. This surface is logarithmically rough in 2- 
d, suggesting that the scaled two-point correlation func- 
tions for the sliding particle and depth problems should 
show an even sharper cusp than for the KPZ surface stud- 
ied in this paper. 

Finally, we remark that it would be interesting to ex- 
plore the effects of removing the hard core interaction 
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between particles. A study of noninteracting particles 
sliding on a KPZ surface in one dimension indicates a 
power law decay of the two-point correlation function 
. A study of other characteristics of the resulting state 
should prove quite revealing. 
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